Hopping and clustering of oxygen vacancies in SrTiOs by anelastic relaxation 
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The complex elastic compliance sn (w,T) of SrTi03_5 has been measured as a function of the O 
deficiency S < 0.01. The two main relaxation peaks in the absorption are identified with hopping 
of isolated O vacancies over a barrier of 0.60 eV and reorientation of pairs of vacancies involving 
a barrier of 1 eV. The pair binding energy is ~ 0.2 eV and indications for additional clustering, 
possibly into chains, is found already at 5 ~ 0.004. The anistropic component of the elastic dipole 
of an O vacancy is AA = 0.026. 



Diffusion and aggregation of oxygen vacancies (Vo) 
in perovskites are still poorly understood, although they 
have implications in several types of applications, like 
in solid state electrolytes for fuel cells or fatigue in 
ferroelectrics. One of the most studied perovskites is 
SrTi03_,5, which has been heavily doped with Vo since 
the 1960s e.g. for inducing superconductivity, but only 
recent experiments fi*2] raised questions on the assump- 
tion that the Vo introduced by high temperature reduc- 
tion are uniformly distributed over the bulk. In addition, 
there is a large spread of the estimates of the activa- 
tion energy for their diffusion J 3 -] , although the prevalent 
opinion is that in most titanate perovskites the barrier 
for Vo hopping is ~ 1 eV p™]. An O vacancy in a 
perovskite lattice has tetragonal symmetry and there- 
fore an associated tetragonal distortion, or elastic dipole, 
which reorients by 90° after each jump causing anelas- 
tic relaxation in addition, being a charged defect, a 
hopping Vo causes fluctuations of the electric field that, 
thanks to the quadrupolar interaction with the surround- 
ing nuclei, may be detected by nuclear magnetic relax- 
ation (NMR). Indeed, two elaxation processes are found 
by NMR in SrTiOs, with activation energies of 0.62 eV 
and 1.1 eV and attributed to hopping of Vo which are 
free or trapped by Fe, respectively [£]. On the other hand, 
SrTiOs reduced in H2 presents a rich anelastic relaxation 
spectrum, including two processes with 0.6 and ~ 1 eV, 
tentatively attributed to H and Vo hopping [2]. Here 
anelastic relaxation measurements on a SrTiOs crystal 
reduced in CO are presented, providing evidence that 
the barrier for hopping of isolated vacancies is indeed 
0.60 eV, as also estimated in recent calculations [2], but 
aggregation into pairs and possibly chains occurs already 
at low 6, so explaining the higher activation energy for 
diffusion generally found. 

The sample was cut from a SrTiOs wafer from M.T.I. 
Corporation as a bar 26.15 x 3.4 x 0.5 mm 3 with the edges 
parallel to the (100) directions. The sample was covered 
with silver paint on a face, suspended on thin thermocou- 
ple wires at the nodal lines and electrostatically excited 
on its first and fifth flexural modes at 5.5 and 74 kHz; 
the electrode-sample capacitance was part of a resonant 
circuit, whose frequency (~ 10 MHz) was modulated by 
the sample vibration, so allowing the vibration to be de- 
tected. The real part of the «u compliance could be 



measured from the fundamental resonance frequency / 
as s' n = p- 1 (0.973fl 2 /hy 2 = 3.64 x 10~ 13 cm 2 /dyn at 
room temperature, where I, h and p are sample length, 
thickness and density [ 7 ]. The elastic energy loss coeffi- 
cient Q^ 1 — s'd/s'u was measured from the free decay 
of the sample vibration or from the width of the reso- 
nance curve. The jumps of Vo cause a reorientation by 
90° of the direction of the nearest neighbor Ti atoms and 
therefore also of the associated tetragonal elastic dipole A 
(actually a quadrupole), having only two independent di- 
agonal elements Ai and A2; the resulting relaxation Ss'd 
of the compliance causes a peak in Q^ 1 (u),T) f£] 
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where AA = Aj — A2, c is the concentration of the relaxing 
defect, vq the molecular volume, to = 2nf, t = Toe w ^ kBT 
is the relaxation rate over a barrier W, the parameter a < 
1 reproduces a possible broadening and the maximum 
occurs at the temperature where lot = 1. 

For the reduction/oxidation treatments the sample was 
inserted in a flattened cylindrical holder of Pt that was 
exposed to a flux of 1 bar of O2 or 0.9 Ar + 0.1 CO 
and heated by induction within a quartz tube cooled 
with water. Each reduction was preceded by oxygena- 
tion for 1.5 h at 950 °C that restored the white translu- 
cent aspect. The O deficiency was deduced from the 
mass change and from the temperature To of the cubic- 
to-tetragonal transformation appearing as a step in the 
elastic compliance near 105 K [ 6 ]. The highest deficien- 
cies after reducing treatments of 3 h at 1100 and 1150 °C 
were estimated from the mass loss as 0.0066 ± 0.001 and 
0.0070 ± 0.001 and resulted in T = 92.4 and 91.7 K re- 
spectively. The assumption of a linear To (d) combined 
with T (0) = 105.2 and 106.9 K before and after the 
cycle of experiments yield T (6) — (106 — 2050 5) K; 
the reported doping 8 is deduced from this relationship. 
The shift of To (0) suggests that the 14 high temperature 
treatments for reduction, homogenization and reoxida- 
tion, each preceded by mechanical removal of the silver 
paste, may have introduced some damage in the sample 
(only part of the results is reported here). 

Figure[T]presents the Q _1 (T) curves in the as-received 
state (curve 0), and after various reduction treatments: 
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FIG. 1: Elastic energy loss coefficient of SrTi03_a measured 

at 5.4 kHz at various concentrations of Vo- The sequence is 
according to S and not chronological. 



S = 0.00088 after 1 h at 950 °C (curve 1), S = 0.0018 
after 1 h at 1000 °C (curve 2), 5 = 0.0041 after 3 h at 
1085 °C (curve 3), 5 = 0.0066 after 3 h at 1100 °C (curve 
4), S = 0.0070 after 3 h at 1150 °C (curve 5). After 
reduction, the sample was cooled to room temperature 
in 2 min (curves 3,4 and 5) or homogenized for 1 h in 
the same reducing atmosphere at 800 °C (curves 1, 2, 
3'); at such low temperature no further O loss occurs. 
The Q _1 (T) curves confirm our previous measurements 
on ceramic samples [ 6 ], with 6 peaks here labeled P1-P6 
starting from high temperature. 

All of them are shifted to higher temperature when 
measured at higher frequency, indicating that they are 
due to thermally activated relaxation processes; this is 
shown for curves 1, 2 and 5 in Fig. O Peak PI has an 
activation energy of ~ 1 eV and therefore it had been as- 
sociated with the jumps of Vo(we were not aware of Ref. 

while peak P3 had been associated with H hopping, 
since its intensity saturates immediately and the reduc- 
tions were carried out in H2 atmosphere [&]. Here there 
is no reason for H contamination, and therefore such an 
assignment must be excluded. Peak P5 is certainly as- 
sociated with doping, but with a very weak dependence 
for 5 > 0.001; it might be related to polaronic relaxation 
and it will not be considered further; a minor peak P6 
and indications for another peak between P3 and P5 are 
hardly detectable. 

The frequency dependence of the anelastic spectra is 
shown for curves 1, 2 and 5 in Fig. [2j and the changes in 
the intensities of the peaks at higher frequency are fully 
meaningful, since both frequencies were measured during 
the same run. The main features are: i) the intensity of 
P3 is oc 1/T at the lowest S, when PI is undetectable, and 
ii) becomes a strongly increasing function of T at higher 
5; Hi) the intensity of PI decreases faster than 1/T. All 
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FIG. 2: Anelastic spectra corresponding to curves 1, 2 and 5 
of Fig. 1, measured on the 1st (circles) and 5th (lines) modes. 



these observations may be explained if P3 is due to hop- 
ping of isolated Vo, while PI is due to the reorientation of 
pairs of Vo- Indeed, the activation energy of P3, 0.60 eV, 
is in good agreement with the NMR experiment [£] , while 
clustering of the Vo in SrTiC>3 has been proposed 
and recent calculations suggest that 2nd neighbor pairs 
of Vo should be stable [2], The proposed picture is as 
follows: at very low O deficiency, practically all Vo are 
free, with a concentration c/ ~ 5; the intensity of PI is 
oc Cf/T and therefore oc 1/T; already at 5 > 0.001 the 
concentration c p of pairs of Vo is significant and, due 
to thermal dissociation, c p (T) is a decreasing function, 
while Cf (T) becomes increasing, since c/ + 2c p ~ 5. This 
explains both the pronounced decrease of the height of 
PI with T and the fact that the height of P3 reverts to 
an increasing function of T. In view of these considera- 
tions, P5 could not be assigned to hopping of isolated Vo 
because its intensity always decreases with temperature, 
whereas it should increase due to the thermal dissociation 
of Vo pairs, when PI is present. 

The first attempts to fit all the spectra have been done 
taking into account only isolated and paired Vo, whose 
concentrations may be estimated as stationary solutions 
of rate equations for the formation of pairs with binding 
energy E p . In this manner, however, it is impossible to 
fit all the spectra with the same set of E p and elastic 
dipole anisotropies AA; in fact, the actual intensity of PI 
increases with 5 much less than predicted by the model. 
This suggests further clustering of Vo, e.g. formation of 
chains along the (100) directions, which are predicted to 
be stable [2]; only the Vo at the ends would contribute 
to PI, since the internal Vo would be bound stronger 
and require a higher energy for jumping out of the chain, 
hence contributing to anelastic relaxation at higher tem- 
peratures. The quantitative description of the formation 
of Vo chains cannot be worked out in a simple manner, 
and it has been chosen to adopt the grandcanonical for- 
malism with the approximation of dividing the lattice 



FIG. 3: a) possible configurations of two or three Vo (empty 
circles) with the respective energies. The Ti atoms are at the 
cell centres and Sr at the vertices, b) potential energy profile 
for various jumps of Vo; below each minimum is shown the 
corresponding configuration. 



into small partitions, within which the formation of pairs 
and chains can be treated exactly [3£]. Figure^) shows 
the configurations of pairs and triplets of Vo that have 
been considered, besides the isolated Vo with site energy 
Ef — 0. A binding energy E p is attributed to both pairs 
of the type Vo— Ti— Vo and in adjacent cell faces; E nn 
to nearest neighbor pairs and E c to each Vo within a 



chain, while the ends of a chain contribute with E, 



The grandpartition function fii] 
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is written for three adjacent 



cells, so that it is possible to include the formation of 
chains of up to four Vo ■ With three cells it is still easy to 
count the multiplicities m a of the possible configurations 
a with n a vacancies and with energy E a by simple 
inspection, and the calculation of the chemical potential 
fx and of the concentrations c a is simple enough to be 
integrated into the non-linear least square fitting routine. 
By defining 6 = e -E P /k B T ^ bc = e -E c /k B T^ x _ e »/k B T ) it 

results: Z = 1 + 16x+ (176 + 67) x 2 + (66 b c + 88b) x 3 + 
(b 6 2 + 246 2 ) a; 4 , where 16a; is the statistical weight of a 
single Vo in one of the 16 O sites of three cells, 176a; 2 
is the weight of the possible 17 pairs with energy E pi 
and so on. The nearest neighbor configuration nn, with 
energy E nn is the intermediate step for the reorientation 
of a pair, but its statistical weight resulted negligible 
and is omitted. The chemical potential /x must satisfy 
the implicit equation 



s _ 3 c _ 3fc B T din Z 
16 16 9/i 



(2) 



where it is recognized that the maximum possible con- 
centration c = 16 over three cells corresponds to S = 3. 
This implicit equation for x has to be solved numeri- 
cally for each T . The decomposition of c into the various 
c a is obtained by keeping track of the contributions of 
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FIG. 4: Fit of the experimental curves for S = 0.007 measured 
at 5.5 and 73 kHz with E p = 0.183 eV and E c = 0.267 eV; 
also shown are the components PI, P2 and P3 at 5.5 kHz. 



the various terms in Z: Cf = (16a; + 134a; 2 + 886a; 3 ), 

c p = W t 17 + ( 66c + 88 ) x + ( b c + A8b ) x2 ] and Cc = 
3^1(3 + 6 c x). 

In Fig. [SJa) is shown the potential profile of the relax- 
ation processes corresponding to PI, P2 and P3, where 
Wi are the barrier heights and Ei the site energies tak- 
ing as zero the isolated Vo- Peak P3 is fitted with Eq. 
(fTj) with c = Cf, and PI with c = c p ; the interme- 
diate relaxation involving the nn pair should have an 
intensity proportional to c p /[T ^cosh 2 '{E2/2k^T)], valid 
for the low concentration limit [ 10 ], and rate r^ 1 = 
To 2 1 e- w »/ kBT cosh(E 2 /2k B T), where E 2 = 2(Wi - W 2 ) 
is the asymmetry between the two states. 

The fit to the case 8 = 0.007 is shown in Fig. g] 
and similar fits are obtained for the other concentrations 
5 > 0.001 assuming binding energies E p — 0.184 eV and 
E c = 0.26 eV. The parameters for the isolated Vo de- 
duced from P3 are: anisotropic component of the elas- 
tic quadrupole (AA) 3 = 0.026, hopping rate with W3 = 
0.60 ± 0.007 eV and r 03 = (5 ± 1) X 10~ 14 s; for the pair 
reorientation (peak PI) (AA)j = 1.87 x (AA) 3 , slightly 
less than for two independent Vo, and the reorientation 
rate has W x = 0.97±0.04 eV and t 01 = (7 ± 4) x 10" 14 s; 
no broadening is found (a > 0.95) for PI and P3. 
The intermediate relaxation P2 has W 2 ~ 0.86 eV and 
E 2 ~ 0.17 eV, but is much broader and more intense 
then expected from the intermediate configuration of re- 
orienting pairs, since a 2 ~ 0.35 and (AA) 2 ~ 5(AA) 3 . 
It is therefore possible that also other configurations, 
e.g. third neighbor Vo, contribute to P2; recent cal- 
culations indeed indicate that the interaction energy of 
Vo may be quite large also at higher distances [2]. With 
AA as small as 0.026, the distortion around a Vo is al- 
most isotropic (for interstitial O in bec metals is AA ~ 1 
p]), as predicted [i 2 .]. Peak P5 has W 5 = 0.43 eV, r 05 = 
(1±0.5) x 10~ 13 s and a 5 = 1. 

Only the spectrum with 5 ~ 8.8 x 10~ 4 (curve 1 of Fig. 
[T]) cannot be reproduced with the same parameters, but 
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requires 5 = 6.4 x 10~ 4 , slightly smaller than estimated 
from the transition temperature T but within the error, 
and E p ~ 0.1 eV, otherwise it is impossible to obtain 
the 1/T dependence of the intensity of P3. A possible 
explanation is that there is a fraction of Vo that are 
strongly trapped at lattice defects; the influence of such 
defects would be particularly evident in curve 1, since 
it has the smallest S and was also the last measurement 
of the series. Indications of lattice damage accumulated 
during the treatments are the presence in curve 1 of a 
broad peak at temperature higher than that of PI and 
the shift of To mentioned above. 

Let us comment on the time for reaching true ther- 
modynamic equilibrium. In fact, after the sample was 
reduced at the highest values of 5 and cooled to room 
temperature within a couple of minutes, the distribution 
of Vo was not in equilibrium. This is demonstrated by 
the effect of a homogenization treatment of 1 h at 800 °C, 
here shown in curves 3 and 3' of Fig. [TJ all the peaks 
remain unaffected except for PI, which decreases. This 
certainly cannot be explained by O uptake during ho- 
mogenization, due to the absence of O2 in the reducing 
atmosphere, nor can it be explained in terms of further 
loss of O, which occurs at higher temperatures and re- 
sults in an increment of the height of PI. It must be 
concluded that the concentration of Vo pairs decreases 
during homogenization due to the formation of more sta- 
ble O clusters. Within the proposed picture, the sup- 
pression of PI is explainable in terms of lengthening of 
the chains of Vo which reduces the fraction of pairs and 
chain ends. 

A seeming inconsistency exists between the fact that 
peaks PI, P2 and P3 imply hopping rates of isolated 
and aggregated Vo exceeding 10 4 s _1 already at 600 K, 
while the aggregation into longer chains near 1100 K re- 
quires hours. The situation is quite similar to that of the 
ordering of the O atoms in the CuO^ planes of semicon- 
ducting YBa2Cu306+a; [ 13 ], which also are perovskite-like 
layers. In that case, the isolated O atoms in the almost 
empty CuO x planes (x < 0.3) hop over a barrier as low 



as 0.11 eV and promptly form stable pairs or short chains 
whose dissociation energy is ~ 1 eV, but reaching equilib- 
rium with longer O-Cu-0 chains requires very long times. 
The explanation proposed for such a behavior is that the 
saddle point for an 2 ~ ion to join another O 2- ion or 
chain fragment is higher than that for hopping away from 
it, due to the electrostatic repulsion. The analogous effect 
is shown for the Vo in Fig. Os), where the saddle points 
involving pairs are higher than those between free Vo- 
If we assume that all the saddle points involving pairs 
are similar, as depicted in Fig. [3Jd), then they are higher 
than in the empty lattice by W\ — E p — W3 ~ 0.19 eV. 
Notice that the rate for forming and dissociating a pair 
is not probed by any of the observed peaks (possibly by 
P2), and it is possible that the enhancement of the corre- 
sponding saddle point is even more than 0.19 eV. In this 
manner it would be explainable why the rate for con- 
structing and rearranging the chains is much slower than 
that for the pair reorientation. The fact that it is diffi- 
cult to assess when the sample is in real thermodynamic 
equilibrium, together with the approximate treatment of 
the lattice as an ensemble of triplets of cells, put some 
limitation on the confidence of the binding energies E p 
and E c derived from the fits of the anelastic spectra. 

In conclusion, the anelastic spectrum of SrTi03_a ex- 
hibits relaxation peaks that are assigned to i ) hopping of 
isolated O vacancies over a barrier of 0.60 eV, ii) reori- 
entation of pairs of vacancies over a barrier of 0.97 eV 
and Hi) jumps involving the intermediate step for the 
pair reorientation and other configurations. Sizeable for- 
mation of pairs starts at S as low as 0.002, and there 
is evidence of further aggregation, possibly into chains. 
The pair binding energy is estimated as 0.18 eV, while 
the binding energy of additional vacancies in a chain is 
~ 0.26 eV. 
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